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Abstract 


A purely ab initio potential energy surface (PES) was refined with selected 32 S 16 02 HITRAN data. Compared to 
HITRAN, the root-mean-squares error (ctrms) error for all J=0-80 rovibrational energy levels computed on the refined 
PES (denoted Ames-1) is 0.013 cm' 1 . Combined with a CCSD(T)/aug-cc-pV(Q+d)Z dipole moment surface (DMS), 
an infrared (IR) line list (denoted Ames-296K) has been computed at 296K and covers up to 8,000 cm 1 . Compared 
to the HITRAN and CDMS databases, the intensity agreement for most vibrational bands is better than 85-90%. Our 
predictions for 34 S 16 02 band origins, higher energy 32 S 16 02 band origins and missing 32 S 16 02 IR bands have been 
verified by most recent experiments and available HITRAN data. We conclude that the Ames-1 PES is able to 
predict 32/34 S 16 02 band origins below 5500 cm' 1 with 0.01-0.03 cm' 1 uncertainties, and the Ames-296K line list 
provides continuous, reliable and accurate IR simulations. The .^-dependence of both line position and line intensity 
errors is discussed. The line list will greatly facilitate SO 2 IR spectral experimental analysis, as well as elimination 
of SO 2 lines in high-resolution astronomical observations. 
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I. Introduction 


Life as we know it cannot exist without several key chemical elements, one of which is sulfur, and the most 
important sulfur-containing molecule in many space environments is sulfur dioxide, SO 2 [1,2,3,4,5,6,7,8,9, 10 ]. 
Determination from experimental or theoretical investigations of fundamental atomic, molecular, and nuclear 
parameters is important for the analysis of NASA's space data obtained by the Spitzer Space Telescope (SST), the 
Stratospheric Observatory for Infrared Astronomy (SOFIA), and the Herschel Space Observatory (HSO), as well as 
other past and future NASA and ESA space science missions. 

The Herschel Space Observatory (HSO) opened an exciting new window into the universe: continuous 
wavelength coverage from 60-670 pm with tremendous sensitivity without any atmospheric interference. This 
spectral range is the playground of molecules, and the potential for significant new discoveries in astrophysics and 
astrobiology is huge. While HSO has completed its mission, analysis of data from the mission continues - for 
example see Refs. [9,1 1,12] for recent studies that identify molecules based on HSO data. 

Of the instruments on Herschel, HIFI covers the spectral range 480-1250 GHz (240-625 pm or 16-42 cm" 1 ) and 
1410-1910 GHz (157-213 pm or 47-64 cm" 1 ), and can achieve very high resolution (circa 1 MHz). With the 
increasingly higher resolution of the latest missions and telescopes, positive identification of molecules will be 
possible as well as the determination of fine details of the environment in which the species exist. The challenge will 
be to separate out the features of interest from the multitude of interfering transitions. The problem is that while only 
a handful of lines are required for the definitive detection of a molecule, a complete characterization of the spectrum is 
required to subtract out the contribution of a molecule that contributes many, many transitions in order to reveal what 
lies beneath. 

The report from the workshop on laboratory spectroscopy in support of Herschel, SOFIA and ALMA clearly lays 
out the problem [13]. "For observations of the dense regions in which star fonnation occurs, the youngest stars and the 
protoplanetary disk material surrounding them possess spectra containing a few molecular species with extremely 
numerous, relatively strong transitions throughout the submillimeter. These species are likely to present a serious 
challenge, in as much as their emission will occupy a substantial part of the spectral range available, thus impeding the 
study of other important species. Thus, the discovery of new species of astrophysical and astrobiological interest will 
be severely hampered unless the spectral lines from these 'weeds' can be removed from the data." This report goes on 
to identify the four most prominent (Class 1) weeds, and five somewhat less prominent (Class 2) weeds. Sulfur 
dioxide is a Class 2 weed. The report further states that the spectra of the weeds and all their isotopologues are 
required, that quantum number assigmnents must be made for all transitions, and that the weeds are expected to be 
associated with temperatures primarily in the 100-300 K range. 

In this work, we are directly responding to the need of observers by combining the state-of-the-art quantum 
chemical method with reliable high-resolution experimental data. While theoretical calculations have been carried 
out for many years, recent advancements in algorithms, computer hardware, and theoretical formulations have made 
the accurate prediction of molecular spectra closer than ever. Recently we reported that the “best theory plus 
high-resolution experiment” strategy has been successfully extended from H 2 0 [14,15], NH 3 [16,17,18,19,20] to C0 2 
[21,22]. Near 7000 experimentally measured C0 2 energy levels up to J= 117 and 14,320 cm 1 were reproduced with 
Grms = 0.0156 cm" 1 . The 12 C 16 0 2 Ames-296K line list predictions near 9000 cm 1 has been verified by recent 
experiments [23] as more reliable than old HITRAN model. [22] Relative to C0 2 , much less high-resolution 
experimental infrared (IR) rovibrational data is available for S0 2 . HITRAN2008 [24] and HITRAN2012 [25] are 
essentially the same for S0 2 which only have 13 vibrational states and 14 associated IR bands in the 0 - 4200 cm" 1 
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spectral range. The purpose of the present study is to provide a high quality line list plus IR intensities for SO 2 that is 
suitable for analyzing astronomically observed spectral data from Spitzer, Herschel, and SOFIA, and for modeling 
higher energy and high-temperature spectra. Note that this SO 2 study is not limited to the far-IR range of HIFI or 
ALMA, but rather we will cover the range 0-7000 cm" 1 in this initial study with .7=0-80 for most bands. The demand 
for quantum number assignments and data for all isotopologues is tailor made for the type of approach that we have 
previously used for the FLO, NH 3 , and CO 2 molecules. 

The primary purpose for spectral data of the weeds is to remove their influence from measured spectra to reveal 
the "flowers", but once given our results, new possibilities exist. For example, one could use the measured spectra of 
SO 2 , which would normally be thrown away, and use it to obtain more detail about isotopic abundances. It is well 
known that isotopic abundances reflect the origin of species, and the weeds with their bright spectra might make a 
particularly sensitive probe to some very interesting processes. 

While our primary purpose has been to generate SCL line lists that will be useful in analyzing high-resolution 
astronomical spectra, we note that SO 2 is of great interest in planetary atmospheres as well, especially Venus and 
exoplanet atmospheres. In fact, SO 2 is the second major component of the Venus atmosphere and the SOIR instrument 
on the Venus Express mission measures high-resolution spectra in the 2400-4400 cm" 1 region, which includes a CO 2 
window. Both the ExoMol database [26] and Exoplanet Characterization Observatory (EChO, an ESA mission 
candidate) [27] consider generation of a S0 3 line list an important need for the community. Additionally, monitoring 
the SO 2 emission generated by power stations on Earth is also a critical focus for environment protection. 
Unfortunately, HITRAN has very little coverage for SO 2 , and thus it is expected that our line lists will also be useful in 
analyzing Venus Express data from the SOIR instrument, the modeling of exoplanet atmospheres, and potentially for 
use in the monitoring of power stations on Earth. 

In this paper, we present a refined spectroscopically accurate potential energy surface (PES) for 12 S 16 0 2 , denoted 
Ames- 1 , as well as an Ames-296K IR line list computed with an ah initio dipole moment surface (DMS) without any 
refinement. Section II describes the theoretical approach adopted in the present work, which is similar to that used 
previously for NH 3 and CO 2 [17,21], so only an abbreviated description is given in Section II. Several technical 
details are described in section III, followed by results and discussion in section IV. A summary and future work is 
presented in the final section. 

It should be noted that this study serves as a basis for future SCL studies, e.g. a far-IR simulation with even higher 
accuracy for both line positions and line intensities, as well as studies that cover a larger spectral range and include 
higher J values, which will be needed to model high-temperature spectra of SO 2 . 


II. Theory and Algorithm 

At NASA Ames, the “Best Theory + High-resolution Experimental Data” strategy has yielded a rms < 0.05 cm" 1 for 
H 2 0 line positions,[14] and better than 95% accuracy for H 2 0 IR line intensities. [28] The water line lists generated in 
Ref. 14 and Ref.28 greatly contributed to the growth and purification of water spectral databases. [15,29] For 14 NH 3 , 
<j nm = 0.02 cm" 1 for 6000 transitions arising from ground state levels, and a rms = 0.04 cm" 1 for 1790 “hot-band” 
transitions arising from V 2 levels. [17, 18] For 12 C 16 0 2 , cr rms = 0.0156 cm" 1 for 6873 purely expt-based rovibrational 
levels and better than 80% intensity agreement for more than 80% of 231 experimentally measured bands. [21] 

However, the key point goes beyond simply reproducing the included experimental data. Much more importantly, 
our focus has been to provide reliable analysis and predictions for unknown bands and isotopologues, as well as 
identifying incorrect assigmnents, with the same accuracy as for directly observed experimental data. This degree of 
unprecedented accuracy (i.e. 0.01 — 0.02 cm" 1 , or at least better than 0.05 cm" 1 ) is critically important as it is required 


4 


to aid in the reliable analysis of high-resolution observational data (for example, see Ref. 1 9 where new assignments 
for NH 3 were aided by our earlier work). More interestingly, a recent CO 2 experiment has confirmed the validity of 
our IR line list at -9000 cnf 1 where HITRAN2008 IR intensities exhibited errors of two orders of magnitude. [23] 
Currently, purely ab initio theory is not capable of providing cr nm < 0.1 cm ' 1 accuracy for the rovibrational spectra of 
polyatomics, except for the two-electron system H 3 . Therefore, high-quality ab initio PES’s need to be refined with 
high-resolution experimental data. Additionally, the DMS needs to be critically examined in order to achieve 
excellent agreement with measured band intensities. Specifically, the shape of the DMS will be different to that of the 
PES, and thus the grid one uses will necessarily be different. 

Many technical aspects are involved in the procedure. One critical part is to solve the full-dimensional quantum 
rovibrational Hamiltonian equations (within the Bom-Oppenheimer approximation) accurately, and efficiently. Both 
rovibrational energies and wavefunctions have to be well converged before, during, and after the refinement procedure. 
Approximations in solving the nuclear Schrodinger equation can easily lead to unpredictable noise in the results, 
which ultimately degrades the reliability and precision of the line positions predicted. The VTET program used in 
this work solves the nuclear Schrodinger equation with an efficiently contracted basis, and converges the rovibrational 
energies well. It then saves the wavefunctions for use in the refinement step. Interested readers are referred to 
Ref.[30] for more details. It should be further noted that even higher convergence criteria should be adopted for the 
wavefunctions to be used in intensity computations, because our goal is to provide reliable predictions for IR spectra 
beyond current experimental limits. 


III. Technical Details 

III.l Ab initio calculations and choice of the starting purely ab initio PES 

We started from singles and doubles coupled-cluster with perturbative triples, CCSD(T) [31] calculations plus 
various small corrections. In total, 489 geometries were selected after our initial PES studies exploring extremely high 
energy regions and the quality of the coverage of the PES, including one -dimensional potential cuts. On every single 
point, CCSD(T) calculations were carried out with the following correlation consistent basis sets: cc-pV(A+d)Z.[32,33] 
aug-cc-pVAZ,[ 34,35 ] X=T,Q,5; aug-cc-pwCVAZ,[ 36 ] A=T,Q; cc-pVQZ-DK[33] for the scalar relativistic 
correction; [37] and the Martin-Taylor basis for core-correlation effects. [3 8] In the end we have 26 global ab initio 
PESs for SO 2 , plus 26 more PES if the scalar relativistic correction is included. In order to find the best starting PES 
to use in the empirical refinement, vibrational variational configuration interaction (Cl) results were compared to a set 
of 125 low-resolution vibrational band origins up to 6900 cm" 1 . The band origin set is summarized in Ref.[39], but it 
is not really a high/medium resolution estimate. Instead, it is essentially tracked back to the Dunham expansion 
formula and anharmonic constants derived from Shelton, Nielsen and Fletcher's experimental work in 1953 [40], 
which were fitted from 17 observed band centers. However, the current upper wavenumber limit of j 2 S 16 C >2 
high-resolution IR data is below 4300 cm" 1 , and the reliability of this reference set is estimated to be a few cm" 1 . We 
conclude it is acceptable for the selection of our purely ab initio PES. 

As described above, J=0 vibrational calculations were carried out on all PESs and the computed ./=() states were 
compared to the 125 vibrational terms from 0 to 6900 cm" 1 . In the end, the CCSD(T)/cc-pVQZ-DK PES fit was 
selected since it had the smallest average deviation: 5.34 cm " 1 and 0.11% for all 125 states. Other choices (based on 
the CCSD(T) method) include: aug-cc-pV(Q+d)Z (6.30 cm' 1 , 0.13%), aug-cc-pwCVQZ (5.40 cm" 1 , 0.11%), 
cc-pV(Q+d)Z+rel (5.90 cm" 1 , 0.13%), aug-cc-pV(5+d)Z+rel (6.89 cm-1, 0.14%), or aug-cc-pwCVQZ (with core)+rel 
(7.14 cm" 1 , 0.14%). All other PESs had relative errors >0.20%. Similar to 12 C 16 02 ,[ 21 ] 32 S I6 02 becomes our 2 nd 
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molecule where a QZ level PES has exhibited the best error cancellation, while a PES that includes extrapolation to 
the one-particle basis set limit [41] yields larger errors. Using the CCSD(T)/cc-pVQZ-DK PES fit, we estimate that 
the main effects being cancelled out include core-correlation and one-particle basis set deficiencies common to 
second-row atoms,[42] while some smaller effects include extrapolation to the complete basis set limit, inclusion of 
diffuse functions in the basis set, and higher-order correlation effects. 


III.2 Least-squares fit and PES representation 

To obtain a highly-accurate, realistic and global PES, we divide the PES into short range terms and long range 
terms, i.e. E S ho rt and E Long . 

v= V Short + E Long- 

We follow previous work on H 2 0[ 14], NH 3 [16,17] and C0 2 [21] regarding the choice of V\ ong and E short formulas, as 
well as the weight function used in the least-squares fit. 400 of 489 final points below 30,000 cm' 1 are more heavily 
weighted: 


5 = [tanh(-0.00045-3-(it -3000Q) + 1. 002 002 00|22. 00200200 
weigh t=s / max (E , 3 0 0 0 ^ 

where E is the potential in cm' 1 related to the C 2v minimum. 

The long-range part, V Long , is based on a simple Morse-potential model plus a simple bending model: 


= £d„(1 - ) 2 + £ D, 2 (l - ) 4 + e 


-0.2-(Ar, ■ +Ar 2 z ) 


( A el Aa i 


+ A 


e2 


A a, 4 ) 


where 

D el = 233,156 cm' 1 , D e2 = 5,250 cm' 1 , A e j = 56,000 cm' 1 , A e2 = 50,000 cm' 1 
P = 1.152733 A' 1 , A r, = ri-r re f, r re /=1.43108 A, Aa\ = cos«zoso-cos( 11 9.3209°) 
r and a represent the S-0 bond length in A and the ZOSO, respectively. The choice of Z) e i and D e 2 are obviously 
much higher than the S=0 bond dissociation energy, mainly in order to ensure that the final PES is globally positive in 
the dissociation region where the short range terms have already been damped out. 

The Eshort tenns are damped out as any S-0 bond dissociates or the ZOSO is strongly bent. The criterion used 
for selecting these damping parameters is the final a rms error for the most important energy region, that is 0 - 30,000 
cm' 1 . Mathematically, V S h or t is given by: 

n = 1 

—d a m ) 2 —d a m /i X(A^) 4 ~d a m j)-(Aai) 2 —d a m ^(A«i ) 4 

f — e /= l i=l 

J damp 


where 


dampl = 1.5, damp2 = 3.0, damp3 = 0.1, damp4 = 0.3, r in A and a in radians 
P = 1.0 A" 1 , Ar, = r, -1.43 108 A, A a\ = cos«zoso-cos(l 19.3209°) 
and P is the permutation operator that ensures that f > [(Ar/)'(Ar 2 ) / ] is totally invariant with respect to the interchange of 
the two O atoms. We note that the PES extends beyond 20,000 cm" 1 , because it is necessary in order to properly 
describe the rovibrational energy levels up to 10,000 cm 1 . That is, the value of vibrational energy levels near 10,000 
cm" 1 are impacted to a small extent by the shape of the PES at higher energies, hence the need for the PES to extend to 
higher energies. 
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The final PES includes 219 coefficients. For the indices i,j, and k, we have 0 < i,k <8 ,j< i, and i+j+k < 12. 
An additional set of 31 boundary points were included in the least-squares fit (with a 1/100 fitting weight) to ensure 
the final fitted PES is globally positive. When fitting to the original ab initio CCSD(T)/cc-pVQZ-DK energies, the 
average fitting deviation (absolute value) for the 393 ab initio energy points below 30,000 cm' 1 is 0.21 cm' 1 , while the 
<j nm error is 0.31 cm' 1 . 

III.3 The Purely ab intio DMS 

On the refined PES, 3928 geometries were randomly generated in the 0 - 60,000 cm” 1 range. The number of 
points selected in each 1000 cm 1 interval monotonically decreases from 126 in 1000 - 2000 cm" 1 to 18 in 59,000 - 

60.000 cm' 1 . From 0 to 30,000 cm" 1 , there are 2881 points. CCSD(T)/aug-cc-pV(Q+d)Z finite-field dipole 
calculations were run on these points and thus there are 3638 dipoles to be used in a least-squares fit. It is well 
known that electrical properties, such as the dipole moment or polarizabilities, require a better description of the 
diffuse part of the wavefunction, which is why the aug-cc-pV(Q+d)Z basis set was used. The dipole moments were 
computed as the energy derivatives with respect to an external electric field. 

Similar to our CO 2 work [21], the DMS was fit according to a pseudo-point-charge-on-nuclei model in order to 
maintain the permutation symmetry of the two O atoms. The same pennutation invariant strategy was adopted 
previously for the FI 5 02 + dipole moment surface fit. [43] The S atom is fixed at the origin of the Cartesian system and 
the point charge on an O atom is expanded in bond length changes A r\. An, and bond angle displacements Aa\ = 
coscuzoso — cos(l 19.32°). Then the sums of the charge-position vector product on the two O atoms are fitted to ab 
initio dipole values. 

V = q 0 r ri + q Q 2 • ri + q s ■ r s = q m ■ n + q Q2 ■ r 2 

969 

q Q \ = A r 2 (1 + cosa zoso ) 

n = 1 
969 

9 , o 2 = 'y, C ijk A/y Arj (1 + cos oc zoso j 

n = 1 

where A r = r-1.43108 A, 0< i+j+k < 16. In total there are 969 Cyk" coefficients for the 1 / 01/02 expansion formula. 
The fitted dipoles are invariant to the pennutation of the two O atoms. The weighting scheme from the PES fit was 
applied in order to focus on the region below 30,000 cm" 1 , for which the average fitting error and average relative error 
are as small as 1.8 x 10" 5 a.u. and 0.010%, respectively. The corresponding root-mean-square (RMS) errors a rms are 

3.0 x 10" 5 a.u. and 0.082 %. 

This purely ab initio dipole surface has been utilized in rovibrational intensity calculations to generate an initial 
IR line list at 296K, denoted Ames-296K. By comparing with available experimental spectra and databases, we can 
find potential defects of the DMS for future improvements. Given that this is a starting DMS, it is fairly accurate, 
while some improvements may be required, especially for future high energy and high temperature simulations. See 
more details in the Results and Discussion section. 

TIT-4 HITRAN2008/20I2 Bands, Intensities, and Energy levels 

The FIITRAN-2008 and HITRAN-2012 databases [24,25] include transition frequencies and intensities of 14 
bands at 296K. Pure rotational bands V2=0 and V 2 =l were taken from the CDMS[44] database, plus 7 bands arising 
from the ground state (GS) and 5 “hot” bands arising from V 2 =l. In the 13 related vibrational states, only the GS and 
V 2 =l pure rotational transitions have the highest J = 99. From ,1=0 to .7=99, there are more than 13,000 rovibrational 
energy levels for 32 S I6 02, but not all of them are appropriate for our analysis. We have to remember that only a small 
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amount of data in CDMS and HITRAN were experimentally measured, while the remaining data were derived from 
experimental models. In our previous l2 C l6 02 Ames-1 PES study, we have shown that some model-based HITRAN 
data contain errors large enough to degrade significantly the empirical refinement. [21] 

First we should note that the two purely rotational bands in HITRAN were completely taken from the CDMS 
database. Here, we first use the CDMS purely rotational energy levels to replace the lower state energies of all other 
HITRAN transitions. Then the HITRAN2008 database for 32 S 16 C >2 was transformed and filtered by a 
multi -determination width requirement (A(E_ max - E min) <0.01 cm' 1 ) to get a “reliable” set of rovibrational energy 
levels. In addition, the ierr parameter in HITRAN for line position uncertainty is adopted for further screening. 
Only those levels derived from ierr>4 transitions are kept. This means the line position uncertainties are smaller than 
lxlO' 3 cm' 1 , which is good enough to compare with most rovibrational energy levels and transitions reported in this 
work. 

In the present work, our calculations and analysis go up to J = 80. This is because it is only the first step for the 
high-temperature high-energy SO 2 IR line list project. J=80 should be good enough for room temperature IR line list 
simulations at 296K, as well as high enough to find potential defects in either the refinement or the line list. For 
,7=0-80, there are 12,441 comparable HITRAN levels. 

The high-resolution data we used in the SO 2 PES refinement were taken from HITRAN, at several discrete J 
values. In order to check the accuracy of the refined PES and line list, our predictions (including both interpolations 
and extrapolations) are compared to additional experimental data beyond HITRAN. After 2008 , the Ulenikov group 
have published five high-resolution experimental papers covering 10 vibrational states: 3 vi,[ 45 ] 3V1+V3 and 
Vi+ 3 v 3 ,[ 46 ] 2V1+V2+V3 and V2+3v3,[47] 3v2 and Vi+ 2 v 2 ,[ 48 ] Vi+3v2, 2vi+V2, and V2+2v3.[49]. Different from 
HITRAN, only transmittance were reported and no absolute line intensities are available. Both published band 
origins and spectra have been compared side-by-side with our predictions and are reported in a later section. 

III.5 Variational rovibrational and IR intensity calculations by VTET 

Rovibrational energy levels, wavefunctions, and transition intensities were computed using the procedures 
described in Ref.[14]. We used hyperspherical Radau coordinates, with the S atom as the central atom, and we used 
the Ilk coupling scheme with maximum bending quantum number of 300. The angular matrix elements of the PES 
were determined analytically after making a 181 term Legendre expansion of the PES. This expansion was carried out 
numerically by Gaussian quadrature, with an average deviation between the re-expansion and original values being 
less than 10' 12 a.u. 

We optimized contracted basis functions for each JPS block. The cut off for solving the one-dimensional 
stretching Schrodinger equations were 0.187 Hartree, i.e. 41042 cm' 1 (1 Hartree = 219,474.6 cm' 1 ) with error criterion 
10' 6 for determining the number of optimized quadrature points. [50] For each K , we kept all bending functions with 
energies below 0.15 Hartree (32921 cm" 1 ), and then we coupled these contracted functions, and kept all with energies 
below 0.12 Hartree (26337 cm' 1 ). The rotation-bending functions were made symmetry eigenfunctions as described 
previously. [14] The determination of the contracted stretching and bending functions was carried out iteratively until a 
self consistent value for the lowest energy for that particular JPS block. The excited functions were virtual levels of 
the ground state, i.e. all the other eigenfunctions of the Hamiltonian for the lowest energy (ground) state. 

Finally we coupled the bending and stretching functions, including all functions whose sum of energies was no 
greater than 0.187 Hartree (41042 cm' 1 ), and extracted all roots with energies up to 0.06 Hartree (13168 cm" 1 ). Full 
pennutation-inversion symmetry was exploited in this final diagonalization step. 


In the intensity calculation, we used 24 point optimized quadrature for the stretches, and carried out the angular 
integrals analytically after making a 60 tenn associated Legendre expansion of the DMS. The expansion coefficients 
were detennined numerically using a 72 point Gauss-Legendre quadrature. 

III.7 Refinement 

Similar to previous work, [17, 21] we allow 21 coefficients (up to the quartic level) of the short-range potential 
expansion to vary. Initially, variations were limited to within a +20% range, the limit we set in our NH 3 /CO 2 
refinements. However, this did not work out very well. The limitations were finally removed and all 21 coefficients 
were allowed to vary freely. This led us to conclude that the Ames-0 PES of 32 S 16 02 is more difficult to refine. 

In those earliest trials, one refinement was carried out using 7=0/4/ 10/20/70 Hamiltonian matrices and related 
“reliable” HITRAN energy levels. Ten Hamiltonian matrices and 30129 eigenvalues were involved. The total 
weighted Orms error was reduced from 3.17 cm' 1 to 0.010 cm' 1 . Compared to the Ames-0 PES, nine coefficients 
were found changed by more than 20%: 521% for A Ra 3 A Rb , -279% for A Ra Ae 3 , -287% for A Ra 2 A 0 2 , 177% for A Ra A Rb , 
121% for A 0 , 120% for A 0 2 , 61% for A R , 42% for A Ra 2 A Rb A 0 , and 41% for A Ra 3 A 0 . Two coefficients for A Ra A 0 3 and 
A Ra “A 0 had their signs changed. Note these relative changes are on short-range terms and they can be large. 
Usually they should be molecule specific and depend on the initial PES. 

Because we were satisfied with this refinement, we computed energies of the fitted grid points on this primitive 
refinement, supplied these to the least-square fitting program together with the 3 1 pre-detennined boundary points. 
The newly fitted PES was used to re-run the 7=0/4/20/50/70 calculations. However, discrepancies of -0.10 cm-1 
were found between some refined energies (based on diagonzalized truncated Hamiltonian matrices) and the new 
variationally computed Cl energies (based on the true exact quantum Hamiltonian within the Bom-Oppenheimer 
approximation frame). We had used 1.4398 A as the reference bond length for R s . 0 in the initial fitting, which is 
0.6% longer than the final value 1.4310 A. This 0.0087 A difference is too large, and could probably account for the 
discrepancies we found here. In order to achieve better consistency, we carried out a 2 nd -round of refinement using 
the new Hamiltonian matrices computed on the primitive refinement. Note that we run the re-fit with its exact 
minimum. 

In this 2 nd refinement based on the re-fit, 5/43/183/158/57 l 2 S l 6 0 2 levels were incorporated with 
2. 5/1. 0/1. 5/2. 0/3.0 weights for 7=0/4/20/50/70, respectively. Ten truncated Hamiltonian matrices (and 38467 
eigenvalues) are computed at 7=0/4/20/50/70. Accordingly, 80/73/49/47/43/30/22/74/28 levels were taken from 
GS/V 3 / 2 V 3 /V 2 /V 2 +V 3 / 2 V 2 /V 1 /V 1 +V 2 states, respectively. An additional 19 high-resolution band origins were taken from 
Ref.45, including 5 repeated values for 2 v 3 , V 2 , 2vi, Vi and V 1 +V 2 . In total, 465 reference values were used. Note 
that lower weights (0.01) were applied on 4/4/2 levels for 7=20/50/70, though this does not mean the 10 reference 
values were wrong or unreliable. The apparent discrepancies were caused by a twisted order of roots. When the 
root order is corrected in the final refinement, the corresponding energies agree with HITRAN values to -0.01 - +0.02 
cm' 1 . 

In the end, the weighted a RM s error was reduced from 0.13 cm' 1 to 0.024 cm' 1 , and the unweighted CTrms was 
reduced from 0.070 cm' 1 to 0.010 cm' 1 . It should be noted that the whole short-range expansion coefficients have 
changed during the re-fit, including higher order terms, so we do not compare to original ab initio values. Because 
the new base PES was fitted to its own minimum, both A r, and A aj gradient coefficients are less than 2E-9 while their 
refined values are 1-2E-5, which are still relatively small. The largest changes are -218% for A Ra A 0 , -52% for 

22 2 3 

A Ra “A 0 , 123% for A Ra “A Rb A 0 , -2743% for A Ra A Rb . Such large relative changes mainly result from their small 
magnitudes, e.g. -0.0192 (A Ra A 0 ), +0.0121(A Ra A 0 ") and -O.53(A Ra ‘A Rb A 0 ). The coefficient for A Ra A Rb varies from 
-0.0041 (Ames-0) -> -0.025 (first refined) -0.0017 (refitted) +0.046 (finally refined). All tests and calculations 
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up to now have not found any unreliability or mistakes that could be associated with these changes of the short-range 
PES terms. 


IV. Results and Discussion 

IV.l Properties of Ames-1: Band origins 

As mentioned above, 8 j 2 S 16 02 vibrational states were included in the refinement. The refined Ames-1 PES 
predicts band origin Go term values with -0.02 cm* 1 uncertainty, for both interpolated bands and extrapolated bands 
beyond 4200 cm' 1 , the EIITRAN limit. See Fig.l, which includes everything experimental we have found that is 
high-resolution. It clearly shows that the accuracy achieved for l4 S 1 ''O 2 is essentially the same as that for 32 S 16 02- 
Similar accuracy for isotopic substitution was found for CO 2 where the 13 C 16 02 Go terms basically have the same 
accuracy as those of 12 C 16 0 2 . 

Full lists of Ames-1 J= 0 vibrational band origins are given in Table 1 for both 32 S 16 C >2 and 34 S 16 02- 
High-resolution experimental data or experimental-model based values were taken mainly from 

Refs. [45, 46, 47, 48, 49, 5 1,52] Note the G 0 term errors for certain bands that are included only for modeling purposes 
(i.e. not actually measured) can be as large as 0.1 - 1.0 cm* 1 , e.g. 2vi+2v2, 4 V 2 +V 3 , vi+4v2. Low-resolution 
“experimental” values (italic in parentheses in Table 1) were taken from Ref.[39] and were computed from Dunham 
expansion formula parameters taken from Ref. [40]. 

The largest deviations for j 4 S 16 02 in Table 1 are the A(Ames-l - Expt) = 0.0222 cm' 1 for vi+v 3 and 0.0487 cm' 1 for 
2v!+V3. The corresponding 32 S 16 C >2 errors are 0.0206 cm' 1 and 0.0255 cm" 1 . For these relatively low-lying vibrational 
states, the relatively bigger discrepancy on 2vi+v3 between 32 S and 34 S isotopologues is suspicious. This assertion is based 
on our experience with CO 2 , NH 3 and the other SO 2 band origins. If the experimentally modeled Go for 2vi+v 3 band 
(3598.7739 cm" 1 ) is verified by future experiments, we may need to investigate further the source of such isotopologue 
dependent discrepancies for the Ames-1 PES. 

[Fig.l Accuracy of Ames-1 PES: 32 S 16 02 and 34 S 16 02 band origins.] 

[Table 1. 32/34 S 16 0 2 band origins up to 5165 cm 4 , computed on Ames-1 PES and compared with available 
experimental data. Note: new 32 S 16 (>2 band origins in Refs.45-49 and 34 S 16 0 2 band origins in Refs.51-52 were 
not included in the Ames-1 empirical refinement.] 

The vi+3v 2 band in 32 S 16 C >2 is unique because its K d = 1 1 levels and transitions are in strong resonance with the 
corresponding 2v 3 K a = 1 1 values. No other K d levels have ever been assigned for this band, which means its 
“experimental” band origin was actually derived from modeling, not directly observed. Extrapolation results from 
K d =\ I to J/K. d =0 largely depend on the specific Hamiltonian formula and spectral lines included. This may explain 
why the vi+3v 2 band origins derived from two high-resolution experimental models differ by nearly 0.10 cm 1 . It 
easily explains the relatively larger uncertainty associated with this transition. Therefore, it is not really a surprise that 
A(Ames-l - Expt) for vi+3v 2 was the largest error for 3 "S 16 C >2 in Table 1, i.e. 0.0511 cm' 1 . The reference value 
adopted in our Ames-1 PES refinement (in 2011), 2693.6348 cm* 1 , was taken from Ref.45, which in 2012 was updated 
to 2693.7279 cm' 1 . [47] If we re-run the empirical refinement with 2693.7279 cm" 1 as the corrected reference, the 
new vi+3v 2 band origin on refined surface comes out at 2693.71 cm* 1 , i.e. now the A is reduced to less than 0.02 cm 1 . 
These results strongly suggest the 2012 value (2693.7279 cm" 1 ) is probably more reliable than the previous one. On 
the other hand, although the refined Ames-1 PES was slightly perturbed by the misleading reference value, its 
self-correcting algorithm still yields a vi+3v 2 band origin toward the right direction. This is a further evidence of the 
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robustness of our refinement procedure. Note that changing that specific reference value used in the refinement stage 
does not have any observable effects on the other HITRAN states because no other vi+3v 2 J/K a > 0 levels were included 
in the fit and no other bands in HITRAN have higher quanta for both vi and v 2 . 

J= 0 levels in Table 1 are cut off at 5165 cm" 1 , which is the highest experimentally related energy. We fill in all 
missing band origins and assigmnents for both 32 S 16 0 2 and j 4 S 16 0 2 so future experimental analyses can benefit from 
our reliable predictions. Average uncertainties for the predictions should be -0.02 cm" 1 or less. Due to the high 
density of states, a full list of both 32 S 16 0 2 and 34 S 16 0 2 J= 0 states from the zero-point energy to 7000 cm 1 is provided 
in the supplementary EPAPS material. [53] We estimate the accuracy of those values will slowly degrade from -0.02 
cm' 1 to 0.10-0.50 cm 1 at 7000 cm" 1 , but only high-resolution experiments with definite low K a lines can verify our 
predictions. 

It should be emphasized that we deem the collaboration with experimentalists extremely important for such high 
resolution IR line list project, because close collaborations involving theoretical computations and experimental 
analyses will greatly facilitate and speed up construction of reliable IR line lists. For example, if 3-5 new band 
origins are experimentally determined near 7000 cm" 1 and become available, performing a new refinement including 
these data will immediately guarantee that all the missing band origins between 5165 cm 1 - 7000 cm" 1 may be 
predicted with better accuracy than we currently have. 

IV.2 Properties of the Ames-1 PES: Geometry, dipole moment, and rotational constants. 

Experimental dipole moment values and formula were given by Patel, Margolese and Dykea in 1979, [54] with a 
reported small uncertainty of 2-4x1 0‘ 5 D. They have been widely accepted in spectral analysis and models, e.g. in 
CDMS. At the Ames-1 PES minimum, our fitted CCSD(T)/aug-cc-pV(Q+d)Z dipole is 1.629402 D, which is 0.16% 
larger than the experimental dipole 1.62673 D. Agreement for the equilibrium structure is excellent, 1.431086 A / 
119.319° (Ames-1) vs. experimental value 1.43076(13) A / 119.33(1)° reported by Patel and co-workers in 1979. [55] 
Both the minimum structure and dipole we report here are in better agreement with experiment than the “best 
estimated” values predicted by Martin et al.[42]. 

For the vibrationally averaged structure and dipole, we are still within 1% errors. For the GS, Vi, V2, and V3 
levels (GS = ground state, i.e., the zero-point level), the experimental dipole values in Ref.[54] are 1.63305 D, 1.63322 
D, 1.62614 D, and 1.65246 D, respectively. We did not compute these vibrationally averaged dipoles through a 
variational Cl approach. Instead, we transformed the Ames-1 PES into a quartic force field (QFF; expanded around 
its exact minimum by 0.005 A/rad step lengths) and use 2 nd -order perturbation theory to compute the vibrationally 
averaged structure first. Then calculate the dipoles of these structures on our Ames-0 DMS. All our dipoles 
computed this way are higher, i.e. 5 (Expt. - this work) = -0.00503D, -0.01182D, -0.00332D, -0.00446D for GS, Vi, V2, 
and V3, respectively. This amounts to -0.3%, -0.7%, -0.2%, -0.3% deviations, respectively. This is not an exact 
comparison since the dipole of a vibrationally averaged structure is different from the vibrationally averaged dipole. 

The vibrationally averaged structure for the GS, Vi, V 2 , and v 3 levels are r(SO)o = 1.435172 A, 1.439009 A, 
1.434929 A, 1.439739 A; and a^oso = 119.344°, 119.392°, 119.517°, and 119.176°, respectively. Compared to 
Martin’s work, [42] we find very consistent geometry changes from the equilibrium structure to the GS structure. For 
example, 5 (ra) = (r,a ) 0 ~(r,a) c : 0.00408A and 0.026° (this work) vs. 0.00397A and 0.02?°. [42] 

Our equilibrium rotational constants agree excellently with those reported in Morino et al.[56] AJBJC e = 
2.01756/0.34556/0.29506 cm" 1 (Expt) vs. 2.01735/0.34546/0.29495 cm" 1 (Ames-1). The reported experimental 
uncertainties are less than 2-3 x 1 O’ 5 cm" 1 . Differences between the two sets of values are as small as 1-2 x 10" 4 cm" 1 , or 
3-6 MHz. However, in 1969 Saito reported a slightly larger A e constant 2.01813 cm" 1 , [55] and Martin’s “Best 
estimate” constants [42] were even larger. The results of this study suggest that the values in Morino et al.[56] could 
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be more reliable, though we note that equilibrium rotational constants are not direct observables, and it would be better 
to compare vibrationally averaged rotational constants. 

The vibrationally averaged rotational constants A/B/C for the 32 S 16 0 2 ground state in CDMS are 60788.550 / 
10318.022 / 8799.754 MHz. [60] Our Ao/B 0 /C 0 via 2 nd -order perturbation theory are 60749.761 / 10315.501 / 
8797.185 MHz, or 2.0263939 / 0.3440881 / 0.2934425 cm’ 1 . The differences, computed as (this work — CDMS), are 
-38.8 MHz (0.06% A 0 ), -2.5 MHz (0.025% B 0 ) and -2.5 MHz (0.03% Co), respectively. The ground state A 0 /B 0 /Coof 
j4 S i 6 0 2 reported in Lafferty et al.[51] are 1.967734 / 0.344188 / 0.292246 cm’ 1 . The corresponding values generated 
in this work via perturbation theory are 1.9667886 / 0.3441028 / 0.2921615 cm’ 1 , which are less by 0.00098 / 
0.000085 / 0.000085 cm' 1 or 0.05% / 0.025% / 0.03%, respectively. 

Next we compare the vibrationally averaged rotational constants for the three fundamental states of 3 “S 6 0 2 . For 
v 2 =l, the CDMS A/B/C constants are 61954.816 / 10320.397 / 8783.856 MHz, and the Ames QFF 2 nd -order 
perturbation theory constants are 61873.402 / 10318.615 / 8781.711 MHz. The differences for B/C are similar to those 
for v 2 =0 (GS), while the difference for A is almost doubled. The differences between the A/B/C constants for vi=l 
and the ground state are 0.001185 cm' 1 , -0.001676 cm' 1 , and -0.001414 cm" 1 , respectively, which are in nice agreement 
with those established in 1993 by Flaud et al [57]: +0.001082 cm' 1 , -0.001662 cm" 1 and -0.001412 cm' 1 . The 
differences between V 3 =l and the ground state AtBtC are -0.020415 cm' 1 , -0.001169 cm 1 and -0.001082 cm' 1 , 
respectively. These also agree well with the Flaud et al [57] values: -0.020710 cm' 1 , -0.001166 cm' 1 and -0.001097 
cm" 1 . These comparisons show that the Ames-1 PES is performing very well. 

Of course, another approach to obtain the vibrationally averaged AIBIC for all vibrational states for the Ames- 1 
PES would be by fitting to the related transitions from the Ames-296K IR line list, though we have not done that here. 

IV.3 Properties of Ames-1 PES: rovibrational energy levels 

The comparison for rovibrational energy levels is based on matching every HITRAN or new experimental 
(denoted NEW-EXP) level to the closest level in the appropriate JPS symmetry block computed on the Ames-1 PES. 
This method works well and no dual matches have been found. As long as the existing assignment can satisfy the 
JPS symmetry, the specific quantum numbers vi/v 2 /v 3 and KJK C are not important. 

Originally there are 13 vibrational states (including the GS) of j 2 S 16 0 2 in HITRAN. But vj+v 3 , vi+v 2 +v 3 and 3v 3 
rovibrational levels are excluded from our comparisons due to the ierr and multi-determination limits, see Sec.IlI-4. For 
new experimental data, the maximum multi-detennination width (uncertainty) is 0.0044 cm" 1 , so all 3959 levels of 7 new 
vibrational states have been included. See results in Table 2. For 12,441 J=0-80 HITRAN-based levels, the deviation 
varies from -0.057 cm 1 to 0.065 cm" 1 , with 0.0001 cm' 1 (mean) ± 0.0122 cm' 1 (Ict). The total root-mean-square 
deviation ctrms = 0.0122 cm' 1 . 

[Table 2 Statistics for J=0-80 energy levels in HITRAN and New Expt. levels, (in cm' 1 )] 

This overall agreement has led us to believe that the HITRAN data included in Table 2 are probably 
self-consistent. Most of the remaining energy levels excluded here have even higher K a , e.g. for GS AT a >3 5 or for v 3 
Kj> 33. For example, the .7=34-45 and 7/ a =34/35 levels of the v 3 band have much larger discrepancies: A(Ames-l - 

HITRAN) = -0.340 cm 1 for the ^ a =34 level, or A = -0.470 cm' 1 for the ^ a =35 level, with an exception A = 0.284 cm" 1 
at J/AT a =3 8/35. Our tests have concluded this is not caused by any problems in our variational VTET calculations. 
This strongly suggests that the /T, -extrapolation degradation in the existing spectral models needs further investigation. 
See Section V for more details. 

It is interesting to note the HITRAN J>50/K !i =\ 1 levels of the v!+3v 2 band and 2v 3 bands need to exchange their 
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band names to recover the energy consistency for each band. Fig. 2 shows the breaks in the existing HITRAN (left, 
before exchange) and the smooth curves at right (after the exchange). This exchange is already included in Table 2. 
In addition, IR intensity tests for the related HITRAN transitions found similar breaks (not shown here). 

[Fig.2 New Assignment for J > 50 / K a =ll levels of Vi+3v 2 and 2v 3 levels: (a) with current HITRAN band 
names; (b) with exchanged band names.] 

In new experimental levels, only one level is confirmed as an outlier, v 2 +3v 3 Jica,Kc = 33 li32 at 4909.9404 cm' 1 , where 
the A(Ames-l - Expt) = 0.315 cm" 1 . We have consistent -0.029 cm" 1 — 0.035 cm" 1 deviations for all other v 2 +3v 3 AT a =0/l 
or J>33 levels. So this level is excluded from the statistics in Table 2. After excluding this level, deviations for the other 
3958 new experimental levels varies from -0.171 to 0.096 cm" 1 , with -0.0191 cm" 1 (mean) ± 0.0268 cm" 1 (la) and a total 
<7rms = 0.0333 cm 1 with respect to zero. All other relatively larger deviations found with the new experimental data (i.e. 
A close to -0.17 cm" 1 or A close to +0.09 cm" 1 ), are associated with the higher end of K d . It can be reasonably explained by 
defects in either the Ames-1 PES or the experimental models, but further tests and or experiments will be needed to 
distinguish between these possibilities. In the future, we will determine whether the higher- A), related errors can be 
minimized in the development of the Ames-2 PES. See Section VI for the updates coming from a new H e ff model. 

Patterns of energy level errors along with the J , K a and energy have been examined for both the HITRAN levels 
and the new experimental levels. See Fig. 3 for the number of levels and Orms for each J, K a and 250 cm 1 interval. 

[Fig. 3. Number of HITRAN and new experimental levels and deviations at each J, I( a and each 250 cm" 1 
interval. Solid circles and triangles are HITRAN levels, while the empty circles and triangles stand for new 
experiment results.] 

Fig. 3 suggests the following conclusions: (1) from 0 to 4200 cm" 1 , the J/K a dependence of energy deviations may 
be further reduced in a future improved version of Ames-1 PES; (2) when extrapolated beyond 4200 cm" 1 , Ames-1 
performs very well for band origins and low J/K a levels, while larger deviations (Ames-1 - New Expt) mainly occur 
for high K a levels. 


IV.4 Ames-296K IR line list for 32 S 16 0 2 
IV.4.A Basics of the Ames-296K IR linelist 

,7=0-80 calculations for 32 S I6 0 2 only printed out all eigenstates up to 0.06 a.u. (13,168 cm" 1 ) including the 
zero-point energy (ZPE) 1535.6336 cm" 1 . An initial IR line list has been generated at 296K with the purely ab initio 
DMS (Ames-0) and all 7=0-80 rovibrational state wavefunctions. The converged total partition function with zero at 
ZPE is 6336.789 at 296K, and 6487.82 at 300K. The corresponding HITRAN partition sum at 296K is 6340.277. [58] 
The difference is ~0.05% and can be ignored in the intensity comparisons. In total, there are 4,094,986 transitions with 
intensity cutoff set to 10" J ° cm.molecule" 1 at 296K. However, it should be noted that this cut-off is not deep enough for 
IR spectral simulations >5500 cm" 1 . Fig. 4 shows the overall comparison between HITRAN and our Ames-296K IR 
line list, with a Gaussian convolution of la = 1 cm" . The cut offs applied here are 10" , 10" " and 10" 
cm.molecule’ 1 . The number of lines in 10" 32 and 10" 36 are 10,890,241 and 63,981,072. Natural isotope abundance is 
used in the comparison, i.e. the Ames-296K is 100% purely 32 S 16 0 2 while the HITRAN intensity values have been 
scaled back to 100% abundance by dividing them by 0.9457. 
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[Fig. 4. Compare HITRAN2012 and Ames-296K Line lists for 32 S I6 C> 2 . (use 100% isotope abundance)] 

Ensuring the convergence of the 32 S 16 C >2 wavefunctions is the most difficult part in performing the SCT IR line 
list computation. The contracted basis sometimes fails to converge under certain conditions, which results in several 
false peaks in the spectra. Occasionally, computational defects affect the energies of lower states, but some do not. 
Experimental/HITRAN data are available only in a very limited energy region, and at limited K a values. Hence, if a 
defect did not affect those HITRAN/new-experiment energy levels, we would not notice it until a simulated spectra is 
plotted using the already computed line list. We have endeavored to eliminate all possible defects in the Ames-296K 
line list, and present it as the best alternative available for bands that have yet to be measured, but it is difficult to 
guarantee that no noise or fake peaks exist within the millions of lines we have computed. We have found that SO 2 is 
even more difficult than CO 2 , due in large part to the density of states. 

A limitation that results from the J and energy cutoffs adopted for computation of the Ames-296K line list is that 
the reliability of the region above 8000 cm" 1 is difficult to estimate. Based on experience, the region below 5500 
cm" 1 will be the most accurate and reliable. Again, based on experience, the region from 5500 to 7000 cm" 1 will 
likely start off with similar accuracy as below 5500 cm" 1 , but will degrade on going to higher energies. This 
degradation is likely to be band-dependent and not necessarily monotonic, but the line list should still be very useful 
for assigning new high-resolution experimental data, similar to what was done previously with our NH 3 and CO 2 line 
lists [19,22] 

On the other hand, based on the comparisons we have performed, we choose to keep this section simple and basic. 
That is, we briefly report what have found and leave out detailed discussions for a future paper wherein we compute 
an upgraded IR line list. We first compare to HITRAN and CDMS, and then we show that Ames-296K based IR 
simulations agree very well with the recent high-resolution experimental spectra reported in Ulenikov et al. 
[Refs.45-49] Note that the symmetric residual 5(I)% = 50% x (I Ames / I obs - I obs / I Am es) is adopted for the comparisons 
of IR intensities. 

IV.4.B Comparison to IR transitions in HITRAN 

The HITRAN2012 database includes 69,356 lines for 12 S 16 02 with both ./, and Jf < 80. Excluding the CDMS 
purely rotational transitions (12,328 at V2=0 and 9,007 at V2=l), there are 48,020 transitions for the other 11 IR bands. 
Table 3 presents a statistical summary of errors for each band including frequency, sum of intensity, and individual 
intensities. Not all the HITRAN transitions were included due to large experimental uncertainties (ierr < 4, as shown 
in the energy level comparisons), band name disagreement, and high end K a levels. To avoid those discrepancies and 
the quantum number twisting problem we identified for a small portion of computed levels, we choose to match our 
computed lines to HITRAN transitions if the following criteria are met: J, and J\ agree, parity (A a +n v3 ) agree, and both 
E\ and frequencies agree to within 0.2 cm" 1 . This rule is not equivalent to the rules we used for energy level 
comparison, and all transitions related to the K d = 11 Vi+3 v 3 and 2v 3 levels are excluded, as well. The total number of 
lines included in this comparison, then is 68,965. 

[Table 3. Statistical summary for 13 HITRAN bands, 130 ^-000 is excluded. ] 

Table 3 summarize the key statistical analyses of the Ames-296K line list vs. HITRAN2012 for 13 bands. 
A+ct(A) for transition wavenumbers and 5 %(I)+ct( 5%) for intensities are the most characteristic numbers that should 
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be examined. Note that detailed plots of the distribution of frequency and intensity errors are given in the EPAPS 
material.[Ref.53], and only a brief summary is presented here: 

(1) Ames-296K intensities are systematically higher than CDMS intensities of the V2=0/l purely rotational bands. 
See further discussion in the next section; 

(2) For Ai f Ai type transition bands, e.g. 2v 2 G v 2 and 2v3(^~GS), intensity deviations 5(I)% start from zero 
and bifurcate into branches with opposite K a linear dependences; 

(3) For B 2 <--> A! type transition bands, e.g. V3 and Vi+v 2 +V3<-v 2 , 5(I)% does not exhibit a bifurcation but 
monotonically drops along K a with the negative slope proportional to the V3 quanta change. As a result, the 
worst intensity agreement in Table 3 is found for 3v3(GGS). 

(4) The 2v 3 band 5(I)% has a wide negative tail resulting from its resonance with vi+3v 2 . The resonance reaches 
its maximum value at AT a = 1 1 which was excluded, as its relative intensity residual 5(I)% exceeds -400%. 

(5) Based on plots of the detailed data, there is not an explicit correlation between error in the transition 
frequency and error in the intensity. 

(6) Obvious frequency error outliers will require further investigation such as diverged A(vi) at K a > 25, irregular 
A(vi+V3) at^ a >19, and A(vi+v 2 +V3<-v 2 ) at K a =\ 8. 

It should be emphasized that these values do not change by more than 5% in the various dipole surface tests, from the 
CCSD(T)/cc-pV(Q+d)Z dipole on the Ames-1 PES grid to the CCSD(T)/aug-cc-pV(Q+d)Z dipole on a much larger 
grid. Least-squares fitting tests also have no observable impact on the intensity comparison results. Given the level 
of theory and the basis set adopted in the ab initio dipole calculation, it is unlikely the level of electronic structure 
theory used is the major source of IR intensity discrepancies. 

IV.4.C Comparison to CDMS 

The 528 (302) experimentally based v?=0 (v 2 =l) purely rotational transition lines available online at the CDMS 
website [59] are actually a subset of the original CDMS fitted 680 (356) lines. [60] The 528 (302) lines are given as 
both line frequency and uncertainty where experimentally determined. The fitted spectroscopic models are then used 
to extrapolate to the full CDMS database including 14754 (9808) lines. Here we compare to both the purely 
Expt/fltted set of 528 lines and the full CDMS, but excluding those J>80 lines. Our Ames IR line list has been 
re-computed at 300K, however, the partition function difference is as large as 1.2%: 6487.82 (Ames-300K) vs. 
6408.32 (CDMS 300K). The CDMS value includes only two vibrational states, v=0 (GS) and v 2 =l, in the partition 
function computation. 

For v 2 =0, the max K d in the purely CDMS fitted set is 23 (except for two lines at K d =26 and K d = 28). Fig. 5 (a) 
and (b) plot both frequency and intensity deviations for Ames-300K vs. the CDMS fitted set (left) and the complete 
CDMS set (right). The frequency deviation range 5(E) is -0.0155 — 0.0203 cm" 1 , with mean 5(E) ± <5rms(E) = 0.0012 
± 0.0045 cm’ 1 . Ames-300K reproduces the intensity of most transitions to a A(I)% value in the +5% range: the 
median A% = -0.32%, the mean A% = 0.19% and a rms (A%) = 1.77%. Note the partition sum difference at 300K is 
not included. With the partition sum adjustment included, the intensity deviation A(I)% range is -4.5% — 10.7%, and 
the mean ± uncertainty is 1.42% ± 1.77%. When K d rises, the R Q branch 5% goes from near zero to negative, but the 
5% for most other branches slightly increases. (See EPAPS) These statistics agree with the dipole comparison result 
where our dipoles computed at vibrationally averaged structures agree with the experimental value to within 0.2-0. 7%. 
Fig. 6 shows how the intensity deviations depend on the AJ and A K. d quantum numbers. These differences may be 
partly due to the neglect of distortion corrections in the CDMS intensity predictions [61], as well as from inadequacies 
in the Ames- 1 based calculations. 
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[Fig. 5 Ames-300K vs. CDMS comparison: (a) V2=0, the pure expt data used in CDMS fit; (b) V2=0, the 
full CDMS; (c) V 2 =l, the pure expt data used in CDMS fit; (d) V 2 =l, the full CDMS.] 

[Fig. 6. Branched comparison for pure rotational band (V2=0) intensity discrepancy, Ames-300K vs. the 
pure expt data used in CDMS fit.] 

The comparison to the full CDMS set where the highest K d > 40 shows significantly larger errors, especially at higher 
K a . The 5% error pattern still depends on different AJ and A K d branches. However, it is interesting to note that both 
the J and K d dependences increase with increasing A K d values. This suggests that the source behind these systematic 
errors is probably a term associated with A K d . For each A K d , the PIQ branch deviations have similarly negative J 
dependence but the R branch deviations exhibit an opposite (positive) J dependence. The highest 5% for the existing 
CDMS data (>60%) occur for A/=l and AAT a =-3, but the steepest J or K„ dependence is found for AK d =7. For those 
K a > 35 transitions, intensity 8 % variations are still consistent with J< 35 transitions, but the frequency errors quickly 
rise from 0.002 cm " 1 up to 0.10 cm" 1 . It should be noted that the upper limit of reliable K d for the ground state levels 
was estimated to be 35 in Muller and Briinken’s original CDMS paper. [60] 

For V 2 =l, the max K d in the purely CDMS fitted set is 21, and the full CDMS predicts up to K d = 39. All error 
patterns are very similar to V 2 = 0 . Frequency errors range from -0.013 cm " 1 to 0.018 cm" 1 , with mean ± la = 0.0005 ± 
0.0047 cm" 1 . The intensity deviations A% are in the -3.1% - 5.3% range, with mean A ± 1<5 rms(A) = 0.9% ± 1.6%. 
This clearly shows that the quantum no. twisting problem does not affect the patterns at all. Note that the partition 
sum difference is already included in these deviations. Again, for the high end K d > 33, the intensity error 5% is still 
consistent with other transitions while the frequency errors increase to near 0.10 cm" 1 . 

IV.4.D Simulated spectra vs. recent high-resolution experiments. 

As the new experimental spectra reported in Ref.45-49 are transmittance, we only compare the overall spectra 
features, matching the peaks, but not quantitatively compare the absolute intensities. 

Figure 7 compares 3vi in the 3412 - 3414 cm " 1 spectral range: experimental spectra (top black curve) reported in 
Ref.48 vs. an IR simulation using the Ames-296K line list with a Gaussian line width 0.002 cm" 1 . We are satisfied 
with the overall agreement. The simulated spectrum is down shifted by 0.01 cm 1 . All the major peaks can be 
easily matched within 0.01-0.02 cm" 1 . However, there are two features (black squares at 3413.05 and 3413.91 cm" 1 ) 
in which the experimental model did not reproduce features that were labeled as impurities. [Ref.48]. The peak at 
3413.05 cm " 1 is probably due to an impurity, but the feature at 3413.91 cm " 1 is reproduced in our simulation as P 'Q40 
of 3vi with intensity 1.25 x 10" cm.molecule’ at 3413.9136 cm" . This is very interesting, because the 
corresponding P ^Q40 transition frequency based on the experimental model is actually in good agreement, i.e. 
3413.9084 cm" 1 . 

[Fig. 7 Ames-296K-based IR simulation (bottom red line) vs. observed 32 S 16 C >2 3vi P branch spectra (top 
black line, reproduced with permission from J. Mol. Spectrosc. 255, 111 (2009). Copyright 2009 Elsevier Inc.)] 

A second example is the 32 S 16 C >2 V 2 +V 3 spectra in the 1870 — 1872 cm " 1 range (see Fig. 8 ). Experimental spectra 
was reported in Ulenikov et al.[46] Again our simulated spectra reproduce most of the main features with reasonable 
intensity agreement and 0.01-0.02 cm " 1 line position errors. The missing peak at 1870.8 cm " 1 is not 32 S I 6 02 , as both 
the experimental model and our IR line list do not contain the feature. There are many “grass-like” small peaks 
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attributed to hot-band 2 V 2 +V 3 f V 2 transitions, and our simulation reproduces these features. 


[Fig. 8 . Ames-296K-based IR simulation (bottom red line) vs. observed 32 S 16 02 V2+V3 spectra (top black line, 
reported in Ref.46), including the hot-band 2 V 2 +V 3 4- V 2 features at right end.] 


V. New Effective Hamiltonian Analysis 

Very recently (while this paper was being written), Ulenikov and co-workers [62] presented their latest 
high-resolution IR spectroscopy measurements in the 1000-1500 cm ' 1 region together with a new Effective 
Hamiltonian (H e g-) model analysis for the ground state (GS) and three vibrational bands in the region: Vi, V 3 and 2 v 2 - 
The T max (K a max ) assigned to their observed transition lines are 89(37), 109(28) and 54(9), respectively. They 
determined that the GS rotational parameters used in the previous literature and databases (including CDMS and 
HITRAN) could not reliably describe the ground state combination differences (GSCD) for those levels with K d 
greater than about 26. The discrepancy between their measured and previously-modeled GSCD increased from 
1.0-1. 5E-4 cm ' 1 ( AT a <2 5 ) to 6.05E-3 cm ' 1 (K a = 34). This was mainly attributed to the fact that the old H eff model was 
fit from K . d max = 23 highly-accurate sub-mm transitions [as quoted from Ref. 62, while the CDMS website data does 
include two transitions with K d = 26/28]. Therefore, Ulenikov and co-workers have developed a significantly 
improved fit including 78 new experimental GSCD values with J max = 53 and A" a max =34. This new model successfully 
reduced the CTrms for the 78 GSCD values (34>V a >29) from 2.69E-3 cm " 1 to l.lE-4cm" 1 . Further, the g R ms for the 
old 149 sub-mm experimental transitions were also reduced from 28.2 kHz to 22.5 kHz. This conclusively shows that 
the new effective Hamiltonian model is superior and more reliable than the previous one, i.e. the one that 
CDMS/HITRAN have been using and we compared to our Ames-1 PES and Ames-296K line list in the sections above. 
In addition, with the resonances among the three states included, 51 newly fitted parameters can successfully 
reproduce more than 4000 rovibrational energies which were derived from 12,131 experimental transitions of Vi, V3, 
and 2 v 2 bands with ctrms around 0.0001 cm" 1 . Thus, the known information for the Vi and V3 bands is doubled. 

[Fig.9 Differences of the GS/V 1 /V 3 / 2 V 2 energy levels and transition frequencies between Ames (This work), 
New Effective Hamiltonian analysis in 2013 (see text for details), and HITRAN data (which is CDMS based): (a) 
and (b) show the energy level differences, (c) and (d) show the transition wavenumber differences.] 

The GS levels derived from this latest H eff model have been pre-screened for higher reliability by Prof. 
Ulenikov’s group, and one set of 2173 GS levels have been kindly shared with us, with K” ydx =45 and .7 max =100.[63] 
These GS levels are combined with 12,131 measured v\Iv^I2\2 transition wavenumbers to get the rovibrational 
energies of the associated V 1 /V 3 / 2 V 2 levels. The numbers of the derived levels are 1863(vi), 1628(v3), and 1305 ( 2 V 2 ), 
respectively. In HITRAN2012, we find 23,504 GS and V 1 /V 3 / 2 V 2 transitions (both ./ Llp and ./| ow <80, and ierr>2) that 
can match with the new H e ff based levels, i.e. both lower and upper levels of those transitions belong to the updated set. 
This is the data used next for comparison in Fig.9. 

There are three sets of data to compare: HITRAN2012, Ames, and New.H e ff.Model. We compute all three AE 

differences and include them in Fig.9, which are AE (New.H e ff.Model - HITRAN), AE (Ames - HITRAN), and AE 
(Ames — New.H e ff.Model). Plot (a) and (b) are for energy level differences, and plot (c) and (d) are for transition 
wavenumber differences. 

It is obvious that the blue triangle (New H e ff vs. HITRAN) outliers and the black square (Ames vs. HITRAN) 
outliers are almost always close to each other, in all 4 plots (a)-(d). As shown in plot (a) and (b), the discrepancies 
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between Ames and the new H e ff levels are insignificant until E > 3500 cm" 1 and K d > 42. Those Ames-HITRAN large 
discrepancy outliers at 2000 cm" 1 and 2500~3500 cm" 1 do not exist in the Ames-New.H e ff.Model comparison. 

Examining plot (b), the agreement between Ames and the new H eff (red circles) at high K d is much better than that 
between Ames and HITRAN. Actually, AEs at K d = 43 are now reduced from 0.439 cm" 1 (Ames - HITRAN) to 0.040 
cm" 1 (Ames - new H e ff). Furthermore, only 6 /f a =44-45 levels have AE(Ames - new Fleff) >0.10 cm' 1 and AE max = 
0.194 cm’ 1 , as shown at the right end of plot (b). On the other hand, the most negative AE(Ames — new Fl e ff) for 
energy levels appear around -0.05 cm "' with K d = 0/1 and / close to 80. The corresponding AE(Ames - HITRAN) 
are similar for those energy levels. However, the most negative AE(Ames — HITRAN) were -0.06 — 0.14 cm" 1 for 
the GS and Vi at K d = 31/32, /=31-40. Most of these larger deviations have been significantly reduced, e.g. the new 
AEs at K d = 3 1 are -0.009 cm" 1 for the GS and -0.039 - -0.047 cm" 1 for Vi. 

In short, the new H eff model and our work show remarkable similarities in the wide range of // =23-42. This 
clearly proves the reliability and higher internal consistency of our “High-resolution Experimental Data + Best 
Theory” approach. As we have stated before, the most valuable part of our approach is to provide “truly reliable” 
predictions (and alternatives) for those unknown or hard-to-measure / analyze spectra. Without knowing any new 
experimental data, our Ames-1 PES based predictions at //, > 30 were much more reliable than any other pre-existing 
models. 

In plot (d), one can see that at around K a = 30 both the AE(Ames - new H e ff) and the AE(Ames - HITRAN) were 
negative, while AE(new H e ff — HITRAN) are nearly zero. This may suggest the discrepancy source could be mainly 
on our end, probably on the data we included in the Ames-1 PES refinement. This is part of the planned 
improvements in a future Ames-2 PES refinement, which definitely will utilize these new H eff based analysis results. 

Based on the Fig. 9 (a)-(d) plots, we estimate the new H eff model GS levels reliably extrapolate to K a = 43-45 (for 
an energy uncertainty less than 0.01-0.02 cm" 1 ). Our Ames-1 PES based predictions at higher K d region can provide 
useful reference for future laboratory work / model development and can be further refined when laboratory data 
become available. 

Finally, we consider this a perfect example for a long-time standing question about the reliability of extrapolating 
to higher quantum number values for H eff based experimental models versus theoretical or variational computations. 
Our basic conclusion is that H eff models are able to maintain higher (1-3 orders of magnitude) accuracy in the range of 
accurately fitted reliable data, however theoretical line lists yield better consistency and reliability for regions far 
above the H eff data region. For both the old and new H eff models of l2 S 16 0 2 , the problematic point for K a is about 10 
— 15 above their respective K d ax . As the new H eff model has substantially improved the accuracy at high K. d , all other 
experimental bands and models should probably be re-analyzed before the next HITRAN or CDMS update. However, 
it should be noted that the intensity discrepancies between the Ames-296 line list and the full CDMS line lists still 
require further investigation. 


VI. Summary and future work 

In this present study, we compute an ab initio PES for SO 2 , refine it using selected HITRAN rovibrational energy 
levels, and then compute an initial 32 S 16 02 IR line list which should significantly improve the ability for astronomers 
to remove the j 2 S 16 02 “weeds” from astronomically observed spectral data. The initial IR line list (Ames-296K) is 
computed using the refined PES together with an unadjusted purely ab initio CCSD(T)/ aug-cc-pV(Q+d)Z dipole 
moment surface. The Ames-1 PES predictions for both interpolated and extrapolated bands are compared to recent 
experiments and have shown satisfactory accuracy and reliability. When compared to the HITRAN and CDMS 
databases, most errors have systematic patterns and probably can be tracked down to one (or a few) K d related terms in 
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either our computations or especially the experimental models for high K d . In some cases, the linear correlation 
between the K d , J and discrepancies depends on AJ and AK a . However, because they are very systematic, we believe 
the ~0.02 cm' 1 and ~10% deviations for line positions and line intensities are adequate to detect certain outliers or 
defects in existing models, and in fact comparison to a very recent study where a new effective Hamiltonian model has 
been developed [62] supports that assertion for line positions. The energy levels computed using the Ames-1 PES 
and the Ames-296K line list are expected to facilitate future IR analyses for both missing bands below 4300 cm' 1 and 
higher energy bands beyond 6000 cm' 1 , as well as for 34 S 16 0 2 . 

In the future, we plan to make the following improvements: (1) incorporate Ulenikov et al’s [Refs. 45-49, 62] 
results to achieve 0.01-0.02 cm' 1 accuracy for everything below 6000 cm' 1 and also extend the prediction reliability at 
higher energies; (2) compute a 34 S 16 0 2 IR line list;[64] (3) investigate the source of intensity discrepancies when 
compared to the CDMS database; (4) include diagonal Bom-Oppenheimer corrections (DBOC) for energy levels and 
line lists of 17/ls O isotopologues; (5) compute and include the non-adiabatic corrections which might help the high-end 
K. d agreement; (6) further improve the convergence of variational Cl calculations to extend the IR line list beyond 
10,000 cm' 1 , and .7>80, for higher temperature IR simulations. Additionally, we may carry out an independent study 
focusing on the far-infrared region to yield microwave spectra for 32/34 S 16 C> 2 , with better than 3-6 MHz and 5-10% 
accuracy for line positions and intensity, respectively. 

The Ames-1 PES, the ab initio dipole surface, and the Ames-296K line list can be downloaded from 
supplementary EPAPS material, [53] or available upon request to authors. 
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Tables 


Table 1. 32/34 S 16 0 2 band origins up to 5165 cm' 1 , computed on Ames-1 PES and compared with available experimental 
data. Note: new 32 S 16 0 2 band origins in Refs. 45-49 and 34 S I6 0 2 band origins in Refs. 51-52 were not included in the 
Ames- 1 empirical refinement. 


32 s 16 o 2 34 s 16 o 2 

Band/cm ' 1 Expt. &# This work Calc-Expt Expt s This work Calc-Expt 

A’ -symmetry 


ZPE 


1535.6336 

v 2 

517.8726 

517.8708 

2v 2 

1035.1264 

1035.1121 

Vi 

1151.7130 

1151.7115 

3v 2 

1 55 1 ,7553 b 

1551.7447 

Vi+V 2 

1666.3348 

1666.3226 

4v 2 

(2066.9) 

2067.7815 

V!+2v 2 

2180.33 12 b 

2180.3074 

2vi 

2295.8083 

2295.8146 

5v 2 

(2582.3) 

2693.6348 11 

2583.2272 

Vi+3v 2 

2693.7279° 

2693.6859 

2v 3 

2713.3826 

2713.3862 

2V!+V 2 

2807.1881° 

2807.1715 

6v 2 

(3095.6) 

3098.0789 

V]+4v 2 

( 3205.8 ) 

3206.4707 

v 2 +2v 3 

3222.9725° 

3222.9510 

2v!+2v 2 

(3318.1) 

3317.9021 

3vi 

3432.2877 d 

3432.2839 

7v 2 

(3608.9) 

3612.3254 

V!+5v 2 

(3717.4) 

3718.6665 

2v 2 +2v 3 

(3730.9) 

3731.9384 

2v!+3v 2 

(3828.1) 

3828.0259 

V[+2v 3 

(3923.8) 

3837.5893 

3v!+v 2 

(3939.9) 

3940.3920 

8v 2 

(4121.3) 

4125.9477 

V!+6v 2 

(4228.1) 

4230.2707 

3v 2 +2v 3 

(4241.5) 

4240.3554 

2V[+4v 2 

(4342.7) 

4337.5544 

Vi+v 2 +2v 3 

(4344.3) 

4343.7922 

3v!+2v 2 

(4446.9) 

4447.8721 

4v, 

(4560.1) 

4561.1041 

9v 2 

(4632.9) 

4638.9188 

Vi+7v 2 

(4738.1) 

4741.2725 

4v 2 +2v 3 

(4748.3) 

4748.2031 

2vi+5v 2 

(4845.0) 

4846.4909 

v !+2 v 2 +2v 3 

(4848.1) 

4849.4206 

2vi+2v 3 

(4953.8) 

4953.5546 

3vi+3v 2 

(4955.1) 

4954.7439 

4vi+v 2 

(5066. 7) 

5065.9680 

10v 2 

(5143.8) 

5151.2137 


1521.1034 


-0.0018 

513.5187 

513.5337 

+0.0150 

-0.0143 

1026.4555 

1026.4359 

-0.0196 

-0.0015 

1144.4786 

1144.4881 

+0.0095 

-0.0106 

1538.7202 

1538.7272 

+0.0070 

-0.0122 

1654.8290 

1654.8172 

2050.4207 

-0.0118 

-0.0238 


2164.5172 


0.0063 


2281.5072 

2561.5218 


0.0511 

-0.0420 


2673.6078 


0.0036 

2679.8009 

2679.7922 

-0.0087 

-0.0166 


2788.6360 

3072.0285 

3182.1012 


-0.0215 


3185.2447 

3295.1355 


-0.0038 


3411.0307 



3581.9311 

3689.9829 

3690.1352 

3801.0247 

3797.0989 

3914.9634 

4091.2121 

4197.3139 

4194.4079 

4306.3153 

4299.2482 

4418.2651 

4533.0415 

4599.8466 

4704.0195 

4698.1296 

4811.0122 

4800.8153 

4906.3149 

4920.9537 

5033.7814 

5107.8021 


A” symmetry 

v 3 1362.0603 1362.0585 -0.0018 1345.0946 1345.0827 -0.0119 

v 2 +v 3 1875.7975 1875.7913 -0.0062 1854.6104 1854.5901 -0.0203 
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2v 2 +v 3 

2388.9155 

2388.9226 

0.0071 


2363.4927 


vi+v 3 

2499.8700 

2499.8906 

0.0206 

2475.8280 

2475.8502 

+0.0222 

3v 2 +v 3 

(n/a) 

2901.4657 



2871.8039 


Vi+V 2 +V 3 

3010.3178 

3010.3149 

-0.0029 

2982.1186 

2982.1057 

-0.0129 

4v 2 +v 3 

( 3413 . 0 ) 

3414.4264 



3379.5302 


Vi+2v 2 +v 3 

( 3520 . 2 ) 

3520.1389 



3487.7569 


2vi+v 3 

3629.7619 

3629.7874 

0.0255 

3598.7739 

3598.8226 

+0.0487 

5v 2 +v 3 

( 3923 . 8 ) 

3924.8035 



3886.6711 


Vi+3v 2 +v 3 

4029.3903 

4029.3743 

-0.0160 


3992.8156 


3v 3 

4054.0012 

4054.0243 

0.0231 


4004.1621 


2vi+v 2 +v 3 

4 136.9341 s 

4136.9158 

-0.0183 


4101.8359 


6v 2 +v 3 

( 4433 . 9 ) 

4435.5888 



4393.2193 


Vj+4v 2 +v 3 

( 4537 . 6 ) 

4538.0261 



4497.2867 


v 2 +3v 3 

4559.4340 c 

4559.4025 

-0.0315 


4505.5413 


2v!+2v 2 +v 3 

( 4643 . 5 ) 

4643.4436 



4604.2444 


3vi+v 3 

475 1 .7 169 a 

4751.6997 

-0.0172 


4713.9476 


7v 2 +v 3 

( 4943 . 3 ) 

4945.7665 



4899.1607 


Vj+5v 2 +v 3 

( 5045 . 2 ) 

5046.0924 



5001.1685 


2v 2 +3v 3 

( 5064 . 9 ) 

5064.2215 



5006.3565 


2v!+3v 2 +v 3 

( 5149 . 3 ) 

5149.3807 



5106.0575 


Vi+3v 3 

5164.8507* 

5164.8370 

-0.0137 


5108.2558 



&32 S 16 0 2 High-res Expt and model values were taken from a Ref.45 and references therein; b Ref.46; c Ref.47 ; d Ref.48 (2009); c Ref.49; 


$ 34gi6Q 2 j-jjgh_ res p; X pt and modeled values from Ref.45, while the data was traced back to Ref. 5 1 and Ref. 52. 

# Low-res modeled values (italic in parentheses) from Ref.39, which included 31 “observed” states, and 84 states derived from Dunham expansion 


formula and parameters given in Ref.40. 
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Table 2 Statistics for J=0-80 energy levels in HITRAN and New Expt levels, (in cm' 1 ) 


Vi V 2 v 3 

E min 

E max 

Jmin 

Jmax 

r/- min 

js max 
^a 

No. 

^min 

^max 

CTrMS 

HITRAN 











0 

0 

0 

1.912 

3958.558 

1 

80 

0 

35 

2217 

-0.057 

0.025 

0.009 

0 

0 

1 

1362.696 

4025.577 

1 

80 

0 

33 

1958 

-0.050 

0.037 

0.010 

0 

0 

2 

2713.383 

4436.386 

0 

76 

0 

23 

1092 

-0.031 

0.034 

0.010 

0 

1 

0 

517.872 

3720.194 

0 

80 

0 

29 

1942 

-0.034 

0.065 

0.016 

0 

2 

0 

1035.126 

2296.506 

0 

62 

0 

20 

894 

-0.019 

0.032 

0.009 

0 

3 

0 

1553.654 

2237.936 

1 

45 

0 

17 

502 

-0.049 

0.009 

0.010 

1 

0 

0 

1151.713 

3301.518 

0 

80 

0 

32 

1606 

-0.051 

0.023 

0.011 

1 

1 

0 

1666.335 

3080.041 

0 

45 

0 

21 

758 

-0.021 

0.007 

0.007 

0 

1 

1 

1876.432 

3964.383 

1 

70 

0 

25 

1425 

-0.039 

0.065 

0.018 

1 

3 

0 

2955.938 

3789.613 

11 

57 

11 

11 

47 

0.023 

0.052 

0.042 


Total 


1.912 

4436.386 

0 

80 

0 

35 

12441 

-0.057 

0.065 

0.0122 

New Expt 











0 

1 

2 

3224.869 

3976.472 

2 

49 

0 

13 

386 

-0.035 

0.005 

0.021 

0 

1 

3 

4563.210 

5071.748 

3 

35 

0 

10 

182 

-0.050 

-0.014 

0.034 

1 

0 

3 

5165.479 

5944.312 

1 

51 

0 

16 

473 

-0.121 

0.065 

0.037 

2 

1 

0 

2809.081 

3418.350 

2 

43 

0 

17 

438 

-0.034 

0.002 

0.016 

2 

1 

1 

4148.472 

4519.131 

4 

26 

1 

11 

115 

-0.047 

-0.016 

0.029 

3 

0 

0 

3434.172 

4782.933 

1 

66 

0 

24 

885 

-0.172 

0.096 

0.046 

3 

0 

1 

4752.344 

5590.652 

1 

53 

0 

16 

478 

-0.159 

0.042 

0.047 


Total 


2809.081 

5944.312 

1 

66 

0 

24 

3958 

-0.172 

0.096 

0.0333 
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Table 3. Statistical summary for 13 HITRAN bands, 130 <-000 is excluded. 


Band in 

Jf 

Y 

F 

F 

. . Freq 

No. 

^min 

Freq 

Freq 

Freq 

Sum-Int 

Sum-Int 

Int 

Int 

Int 

Int 

HITRAN 


^min 

Mnax 

^max 

^mean 

±a(A) 

HITRAN 

5 

^min 

^max 

Savg 

±o(8) 

000^000 

0-80 

0-42 

0.0 

265.9 

12291 -0.026 

0.109 

0.003 

0.010 

2.208E-18 

2.3% 

-ii.i% 

64.2% 

10.1% 

11.2% 

010^010 

0-80 

0-40 

0.0 

201.9 

9007 -0.032 

0.130 

0.001 

0.009 

1.779E-19 

2.4% 

-5.4% 

34.6% 

5.8% 

6.0% 

010^000 

0-70 

0-26 

436.6 

645.6 

5916 -0.024 

0.037 

0.001 

0.010 

3.707E-18 

-2.1% 

-55.1% 

32.9% 

-3.3% 

18.9% 

020<-010 

0-62 

0-21 

446.4 

622.1 

3727 -0.019 

0.016 

-0.006 

0.006 

5.770E-19 

-2.6% 

-44.8% 

28.1% 

-3.1% 

16.2% 

030 <-020 

1-46 

0-17 

463.1 

598.3 

1532 -0.032 

0.017 

0.006 

0.008 

5.588E-20 

-3.3% 

-35.5% 

20.3% 

-3.7% 

12.6% 

100^000 

0-80 

0-32 

1031.0 

1273.2 

8256 -0.038 

0.056 

-0.002 

0.005 

3.322E-18 

3.0% 

-47.8% 

28.7% 

0.7% 

11.1% 

110<-010 

0-46 

0-22 

1047.9 

1243.8 

4043 -0.026 

0.019 

-0.005 

0.006 

2.506E-19 

3.1% 

-81.1% 

21.5% 

0.3% 

17.9% 

001^000 

0-80 

0-33 

1298.0 

1405.7 

5427 -0.013 

0.025 

0.000 

0.004 

2.573E-17 

2.3% 

-44.8% 

8.4% 

-3.7% 

6.7% 

011<-010 

0-71 

0-25 

1302.1 

1397.0 

3948 -0.018 

0.049 

-0.001 

0.009 

2.022E-18 

3.8% 

-20.9% 

30.0% 

-0.3% 

6.7% 

101^000 

0-80 

0-24 2433.2 2533.2 

4029 -0.084 

0.110 

0.009 

0.022 

5.434E-19 

-6.4% 

-52.8% 

15.3% 

-11.0% 

5.2% 

111<-010 

0-61 

0-21 

2441.2 2521.1 

2733 -0.144 

0.017 

-0.009 

0.014 

4.239E-20 

-5.2% 

-36.9% 

-1.6% 

-8.2% 

4.1% 

002 <-000 

0-76 

0-24* 

2599.1 

2787.9 

4400 -0.042 

0.011 

0.003 

0.008 
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Figure Captions and Figures 


[Fig.l Accuracy of Ames-1 PES: 32 S 16 C >2 and 34 S 16 02 band origins.] 

[Fig. 2 New Assignment for J > 50 / K a =ll levels of vi+3v 2 and 2v 3 levels: (a) with current FIITRAN band 
names; (b) with exchanged band names.] 

[Fig. 3 Number of FIITRAN and new experimental levels and deviations at each J, Ka and each 250 cm ' 1 
interval. Solid circles and triangles are HITRAN levels, while the empty circles and triangles stand for new 
experiment results. (Refs. 45-49)] 

[Fig. 4. Compare HITRAN2012 and Ames-296K Line lists for 32 S 16 0 2 . (use 100% isotope abundance)] 

[Fig. 5 Ames-300K vs. CDMS comparison: (a) v 2 =0, the pure expt data used in CDMS fit; (b) V 2 = 0 , the full 
CDMS; (c) V 2 =l, the pure expt data used in CDMS fit; (d) V 2 =l, the full CDMS.] 

[Fig. 6 . Branched comparison for pure rotational band (v 2 = 0 ) intensity discrepancy, Ames-300K vs. the pure 
expt data used in CDMS fit.] 

[Fig. 7. Ames-296K-based IR simulation (bottom red line) vs. observed 32 S 16 C >2 3vi P branch spectra (top black 
line, reproduced with permission from J. Mol. Spectrosc. 255, 1 1 1 (2009). Copyright 2009 Elsevier Inc.)] 

[Fig. 8 . Ames-296K-based IR simulation (bottom red line) vs. observed 32 S 16 0 2 V 2 +V 3 spectra (top black line, 
reported in Ref.46), including the hot-band 2 v 2 +V 3 V 2 features at right end.] 

[Fig. 9 Differences of the GS/vi/v 3 / 2 v 2 energy levels and transition frequencies between Ames (This work), 
New Effective Hamiltonian analysis in 2013 (see text for details), and HITRAN data (which is CDMS based): (a) and 
(b) show the energy level differences, (c) and (d) show the transition wavenumber differences.] 
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Fig. 1 Accuracy of Ames- 1 PES: 32 S 16 02 and 34 S 16 C >2 band origins. 



Fig. 2 New Assignment for J > 50 / K a =l 1 levels of vi+3v 2 and 2v3 levels: (a) with current FIITRAN band names; (b) 
with exchanged band names. 
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Fig. 3 Number of FIITRAN and new experimental levels and deviations at each J, Ka and each 250 cm' 1 interval. Solid 
circles and triangles are FIITRAN levels, while the empty circles and triangles stand for new experiment results. 
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Fig. 4. Compare FHTRAN2012 and Ames-296K Line lists for 32 S 16 C> 2 . (use 100% isotope abundance) 
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Fig. 5 Ames-300K vs. CDMS comparison: (a) V2=0, the pure expt data used in CDMS fit; (b) V2=0, the full CDMS; 
(c) V2=l, the pure expt data used in CDMS fit; (d) V2=l, the full CDMS. 
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Fig. 6 Branched comparison for pure rotational band (v2=0) intensity discrepancy, Ames-300K vs. the pure expt data 
used in CDMS fit. 
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Fig. 7. Ames-296K-based IR simulation (bottom red line) vs. observed 32 S 16 02 3vi P branch spectra (top black 
line, reproduced with permission from J. Mol. Spectrosc. 255, 1 1 1 (2009). Copyright 2009 Elsevier Inc.)]). 



Fig. 8 . Ames-296K-based IR simulation (bottom red line) vs. observed ’’S , 6 0 2 V 2 V +3 spectra (top black line, 
reported in Ref.46), including the hot-band 2 V 2 V+ 3 V f 2 features at right end. 
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Fig. 9 Differences of the GS/V 1 /V 3 / 2 V 2 energy levels and transition frequencies between Ames (This work), New 
Effective Flamiltonian analysis in 2013 (see text for details), and HITRAN data (which is CDMS based): (a) and (b) 
show the energy level differences, (c) and (d) show the transition wavenumber differences. 
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